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Abstract 

The most general form of Hamiltonian that preserves fermionic 
coherent states stable in time is found in the form of nonstationary 
fermion oscillator. Invariant creation and annihilation operators and 
related Fock states and coherent states are built up for the more gen- 
eral system of nonstationary forced fermion oscillator. 

PACS numbers: 03.65.-w, 03.65. Ca, 03.65.Vf, 05.30.Fk 



1 Introduction 

The time evolution of coherent states (CS) has attracted a great deal of 
attention since the introduction of Glauber's CS of the harmonic oscillator 
[H [2], [3]. Of particular interest has been the determination of the Hamil- 
tonian operator for which an initial coherent state remains coherent under 
time evolution. It is established that this Hamiltonian has the form of the 
nonstationary bosonic forced oscillator Hamiltonian [H El El 0, IB] : 

H,, = ujit)a^a + /(t)at + /*(t)a + (1) 
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where uj{t) and are arbitrary real functions of time t, and f{t) is arbitrary 
complex function. 

Our purpose in the present article is to study the dynamical invariants and 
time evolution of CS for general (one mode) fermionic Hamiltonian and to 
establish the most general form of Hamiltonian which preserves the fermionic 
CS under the time evolution. 

The organization of the article is as follows. We start with a review in 
Sec. II of some main results of time evolution of bosonic forced harmonic 
oscillator. In Sec. Ill we study the temporal stability of fermionic CS and we 
show, by using the fermionic analog of the invariant boson ladder operator 
method [9|, [TOl [11] , that the most general form of Hamiltonian that preserves 
fermionic CS stable in time is in the form of nonstationary fermion oscilla- 
tor. In Sec. IV we treat the more general system of nonstationary forced 
fermion oscillator (FFO), which is shown to be the most general one mode 
fermionic Hamiltonian system. Following the scheme related to the boson 
system ^ we construct the dynamically invariant fermion ladder operators 
and related Lewis- Riensenfeld Hermitian invariant [12]. Using these invari- 
ants, we construct fermionic Fock states and CS for FFO system, which can 
represent (under appropriate initial conditions) the exact time-evolution of 
initial canonical CS. Finally the relation of the invariant ladder operators 
method P [lO] to the Lewis- Riesenf eld method [12j is briefly described on 
the example of FFO. The paper ends with concluding remarks. 



2 Canonical CS and their temporal stability 

The standard boson coherent states (CS) (called also Glauber CS, or canoni- 
cal CS) are defined as the right eigenstates of the boson (photon) annihilation 
operator a [H [21 E] 

a\z) = z\z) 

the eigenvalue z being a complex number. The annihilation and creation 
operators a and satisfy the commutation relations[a, a^] = aa^ — a^a = 1. 
The normalized CS \z) can be constructed in the from of displaced ground 
state |0) [Il[2l[3|, 

\z) = D{z)\0), D(z) = eH-^*'^), (2) 
and their expansion in terms of the number states \n) reads 

I |2 CO 

k)=e-^E^I^) (3) 
n=ovn\ 
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The problem of temporal stability of bosonic CS is solved by Glauber [1] 
and Mehta and Sudarshan [5] (in the case of one mode CS, and for n-mode 
CS - by Mehta et al. [6]). The result is that the most general Hamiltonian 
that preserves an initial CS \z) stable in later time is of the form of the 
nonstationary forced oscillator Hamiltonian Hcs-, eq. ([T]). The Hamiltonian 
(PP) that preserves CS stable is shortly called coherence Hamiltonian. Thus 
the boson coherence Hamiltonian takes the form of a non- stationary forced 
oscillator Hamiltonian. Here "stable" means that the time evolved state 

i—\z]t) = H^s\z]t) . (4) 

remains eigenstate of a possibly with a time-dependent eigenvalue z{t), 

a \z] t) = z(t) \z; t) (5) 

From the latter equation one deduces that, up to a time- dependent phase 
factor exp{iip(t)), the time-evolved CS \z;t) depends on time t through z(t), 
that is 

\Z; t) = e'^W ^ |^(^)) = ^aU(t)-z'(t)a (g^ 

One says that for system with Hamiltonian ([T]) an initial CS remains CS 
all the later time [U [5] (or remains temporally stable). For the Hamiltonian 
system ([T]) the time dependent eigenvalue value z{t) obeys the equation [HE] 

iz = uj{t)z + fit) 

the solution of which takes the explicit form [z = z{0)) 

z{t) = mz + l{t), ^(t)=e-^/o-(*')'^*', (7) 
7(t) = /e'i'o'^("^''V(iVi' 1 e^'-^o-^*')'^*' (8) 



In the particular case of constant u we have 

z{t) = e-'"«* (z-ij e'^^^'Fit^dt 



The forced oscillator system ([T]) admits linear in terms of a and invariant 
boson annihilation operator Ac{t), [A{t) , A'' (t)] = 1, 

A{t) = U{t)aU\t) = (3{t)a + -f{t) = A„ (9) 
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where U{t) is the unitary evolution operator, and 



For any system the time-evolved CS \z; t) are eigenstates of the corresponding 
invariant annihilation operator A{t) with constant eigenvalues 2, A{t) \z\t) = 
z\z;t), and can be represented in the form of invariantly displaced time- 
evolved ground state \0;t) = U(t) |0) , 

\z;t) = D{z,A{t))\0;t) , D{z, A{t)) = e'^'^'^'~'''^^'\ (10) 

If A{t) is invariant then A^{t) also is, and any other combination of them 
is also invariant. In particular A''{t)A{t) and D{z,A{t)) are also invariant 
operators of the forced oscillator ([T]). Invariant operators are very useful, 
since they transform solutions into solutions, as demonstrated in ffTUj) . 

The invariant boson ladder operator (I9l) is a simple particular case of 
linear invariants of general quadratic quantum system, constructed first in P, 
[To] . For the nonstationary quantum oscillator Hermitian quadratic in a and 
invariant was constructed and studied by Lewis and Riesenfeld [T^. Using 
these properties of the invariants it was shown [TT] that a given Hamiltonian 
H preserves the temporal stability of CS \z) if and only if it admits invariant 
of the form Ac = l3{t)a+'y{t). The general form of such Hamiltonian coincides 
with Glauber- Mehta-Sudarshan coherence Hamiltonian ([T]). 

3 Temporal stability of canonical fermion CS 

Fermion coherent states (CS) are defined (see [13l [TH [IHl [16] ) as eigenstates 
of the fermion annihilation operator b, 

& 10= CIO, (11) 

where the eigenvalue is a Grassmannian variable: = 0, CC* + C*C = 0. 
Recall the fermion algebra: 

{b,b^} =bb^ + b^b = l, b^ = b^^ = 0. (12) 

For definiteness eigenstates of fermion ladder operator b should be called 
canonical fermion CS. This is in analogy to the eigenstates of boson annihi- 
laion operator a, which are known as Glauber CS and canonical boson CS as 
well. In terms of the Grasmann eigenvalues ( many of the properties of \() 
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repeat the corresponding ones of the bosonic CS 1^;) [16\. In particular one 
IC) = ^(C)|0) = e-^^*^(|0)- C|l)). (13) 
rfCXlC)(Cl =1, (14) 

where -D(C) = exp{b^( — C*b), |0) is the fermionic vacuum, 6|0) = 0, and 
|1) is the one-fermion state, |1) = |0). The integrations over ( and (* are 
performed according to the Berezin rules [IE] 



I dCdCCC = 1, y dCdCC = J dCdCC = J dCdCl = 0. (15) 

The temporal stability of the canonical fermion CS is defined in analogy 
to the temporal stability of canonical boson CS, namely the evolution of an 
initial \0 is stable if the time-evolved state \C',t) = U(t)\() {U{t) being the 
evolution operator of the system) remains eigenstate of b in all later time, 

b\c:t) = m\c:t). (16) 

It is clear that the time-evolved states \C'-,t) also obey the overcompleteness 
relation (1141) and are eigenstates of the invariant ladder operator B{t) = 
U(t)bW(t). This means that the B{t) and b should commute (we suppose 
that C(t) and ( commute). The general form of a fermionic operator is 
a (complex) linear combination of b, 6^ and 6^6. Such a combination will 
commute with b under certain simple restrictions. Taking into account that 
the invariants B{t) and -B^(t) have to obey the fermion algebra ffT^ we derive 
that [b,B{t)] = if and only if B{t) is proportional to6, B{t) = [3{t)b. Thus 
the fermion coherence Hamiltonian should admit invariant of the form 

Beit) = pm (17) 

where f3{t) may be arbitrary complex function of time. As we have already 
noted at the end of the preceding section, similar form of the ladder operator 
invariant A^., eq. ([9]), is required in the case of boson systems [11]. To obtain 
now the general fermion coherence Hamiltonian Hfcs we apply the defining 
requirement for quantum time- dependent invariants B{t), 

^B{t)-t[Bit),H] = (18) 

to the operator (JT71) . The general form of fermionic (one- mode) Hamiltonian 
is a Hermitian linear combination of b, 6^ and 6^6, 

Hf = Luitpb + fitp + rit)b + git), (19) 
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where uj{t) and g{t) are real functions of time. The substitution of this Hf 
into (fTSjl for Bc{t) produces the two conditions 

P = i(3uj, = Pf. (20) 

These simple conditions are readily solved, f{t) = 0, = exp (^i J* uj{T)dT^ 
leading to Hamiltonian 

Hfcs = uji0b + git), (21) 

which is the most general form of fermion coherence Hamiltonian. If the 
evolution of an initial CS \() is governed by Hjcs, then the time-evolved 
state IC; t) remains eigenstate of b with eigenvalue 

C(t)=r'(t)C = e-^-/"o-M'^X (22) 

The results fl2ip and (1221) are similar in form, but not identical, to those 
for the boson systems ([T]) and ([7]). The fermion coherence Hamiltonian (12T]) 
is of the form of an oscillator with time dependent frequency {nonstationary 
fermion oscillator) , while the boson coherence Hamiltonian ([T]) is of the more 
general form of the nonstationary forced oscillator. In the next section we 
find the exact evolution of fermion CS and fermion number states, governed 
by the nonstationary forced oscillator Hamiltonian using the time-dependent 
integrals of motion method [9l [TOl fT2] . 



4 FFO and invariant ladder operators 

We consider the single nonstationary fermionic forced oscillator (FFO) de- 
scribed by the Hamiltonian (fT9l) . As we have noted this in fact is the most 
general Hamiltonian of single fermion system. The fermion number operator 
is defined as = 6^6. It obey the relation A^^ = A^ and the three operators 
b, 6^ and A^ close under commutation the algebra 

[6, A^] = b, [b\ N] = -b\ [b, b^]=l- 2N, (23) 

The Hilbert space Ti. of the single-fermion system is spanned by the two 
eigenstates {|0) , |1)} of A^: 

6^6 \n) = n\n) , n = 0, 1 

The operators b and b'' allow transitions between the states as 

6|0) = 0, = |0) , 6^ |1) = 0, 6^10) = (24) 
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Let us also note that linear combinations ofb\ b and produce the half-spin 
operators Jj, 

Ji = i(6t + 6), 72 = ^(6^-6), Js = b^-l 

closing the su{2) algebra: [Jk, Ji] = i^kimJm- It is clear that the fermion 
forced oscillator Hamitonian Hjcs belongs to the central extension of su{2) 
(is a linear combination of Ji plus free C-number term). 

It is convenient to use raising and lowering operators J± = Ji ±iJ2 which 
satisfy the following commutation relation: [J+, J„] = 2J3, [J3, J±] = ±J± , 
where = b\ J_ = b. So that in terms of these half spin operators the 
Hamiltonian ( fTOl) takes the form 

Hf{t) = uj{t)J, + /(t) J+ + r (t) J„ + g{t) + ^. (25) 

Our task is the construction of the time- dependent invariants for the 
system (1191) . 0251) . The defining equation of the invariant operator B{t) for 
a quantum system with Hamiltonian H{t) is (1181) . Formal solutions to this 
equation are operators B{t) = U(t)B{0)U'''{t), where U{t) is the evolution 
operator of the system, U(t) = Texp[— z j ^H{t)dt]. In our case of FFO 
(fT9l) . ( l25l) we look for the non-Hermitian invariants of the form of linear 
combination of the SU{2) generators fl^^ . 

Bit) = iy4t)J_ + u+{t)J+ + p;{t).h, 
B\t) = u*_{t)J+ + u\{t)J^ + ul(t).h, 

where ^±{1), v^{t) may be complex functions of the time. Hermitian invari- 
ants then can be easily built up as Hermitian combinations of B and B'^ . In 
particular if 5 is a non-Hermitian invariant the operator 

/ = B^B - \ (27) 

is an Hermitian invariant, the fermion analog of the Lewis- Riesenfeld quadratic 
invariant [12] . 

Substituting (PUI) . into (|T^ . we find the following system of differ- 
ential equations for the parameter functions of B{t) 

z>3 = 2z(z/+r-z/_/), (28) 
z>+ = i{u-if - u+uj), (29) 
i^- = t{y^uj-u^n. (30) 

The solutions of the above linear system of first order equations are uniquely 
determined by the initial conditions i^±(0) = i^o,±5 ^^3(0) = z^o,3- If we want 
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the invariants B{t) and -B^(t) be again a fermion ladder operator, i.e. to obey 
the conditions 

B\t) = 0, {B{t),B\t)} = l, (31) 
we have to take z/o,± and z/0,3 satisfying 

= -4z^o,+'^o Wo- \ + Wo,+ \ = 1 • (32) 

Indeed, for B'^{t) and {B{t),B^{t)} we find 

B^t) = u+u^ + \ul = Xi{u±,U3), 
{Bit), B\t)} = |z._p + p + i|z.3p = A2(z^±, Z.3). 



(33) 



The quantities Ai(z/-|-, z/3), A2(z^±, 2^3) turned out to be two different ^constants 
of motion^ for the system ( l28|) - (l30l) . their time derivatives being vanishing: 



iA, = i(i..p+Kp+ii.3n=o. 

Therefore we can fix the values of these constants as Ai = 0, A2 = 1, i.e. 

^+^- + i'^3 = 0> k+P + iksP = 1, (35) 

and satisfy the conditions flHTl) . If the initial conditions are taken as 

z/_(0) = l, z/+(0) = = z/3(0), (36) 

then B{0) = h. 

Let us first note that in the particular case of the free oscillator, f{t) =0, 
the solution of the above system of equations is readily obtained in the form 

^.±(^) = i.o,±e±^^*"(^)'^ z/3 = i^o,3, (37) 

where i'o^±, i^o,3 are constants. For this solution the expressions Ai(z/-|-, 1/3), 
\2{i^±, i'^) are readily seen to be constant in time as expected: Ai = i/o,-i^o,+ + 
vl^^/A, A2 = |i^o + ? + |i^o,3lV2- Then from ([M]), ([SSD we see that the 
invariant fermion annihilation operator B{t) now takes the form B^oit), 

BUt) = iyo,-e-'''^% + uo,+e'^^'¥ + 2^-uo,-Uo,+ {b^b - i) , (38) 

where ip{t) = uj{T)dT and |z^o,-| + I'^o.+ I = 1; the phases of z/o,± remaining 
arbitrary. 

Consider now in greater detail the nonstationary forced oscillator with 
nonvanishing f{t)\ f{t)j^O. 
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For this system we can express all the three parameter functions i'±(t), 
1^3 (t) in terms of one of them, which has to obey more simple second order 
differential equation. Let for example, express t^sit) and z^-(t) in terms of 
z/+(t) and its derivatives. We have 



1 

2/2 



1^3 



1)4. + (iuj 



f 



(39) 
(40) 



Taking the time derivative of both sides of (HUj) and using eqs. 
we arrive to the third order equation for z^+(t). 



and 



1 

2/2 



2/3 



2/2 





+ 2/2 


2/2 2/4 


■ojf 
«/3 


2/3) ' 


2/3 


r/ 
p 


■"iCof 1 01) 
2/3 2/2 


iuif 
- 2p 


-r f 2/4 



(41) 



Furthermore if we could find a first integral of the equation fBTl) then we 
can express z/_ in terms of z>+, eliminating the second derivative in eq. 
( HOl) . One can check that the following expression of z/+, z>+ and i>+ is a first 
integral of equation fHTj) (that is dX/dt = 0), 



A 



/2 



(42) 



We regard this formula as second order equation for u^, depending on an 
arbitrary constant A. Using this, and supposing that z/+ 7^ 0, we obtain for 
i>_ the more compact expression in terms of and z>_|_. 



(43) 



and we see that z/_ = (A/4 — z/3)/4z/_|_. The first integral A of eq. (HTj) is 
proportional to the constant of motion Ai of system (l28l) - ( |30l) : Ai = A/16. 

Thus the operators B{t), B^{t), eq.(l26l), are invariant for the forced os- 
cillator (|T9|) . fl25l) if i^+(t) is a nonvanishing solution of the second order 
equation with any constant A, and z/3 and z/_ being given by eqs. flHUj) 
and fHHl) respectively. They will obey the fermionic ladder operator condi- 
tions (l3Tl) if Ai = = A and A2 = 1. Instead of fixing A2 = 1 we can redefine 
B{t) — > B{t) I a/A2, i.e. take the invariant fermion annihilation operator of the 
form valid for any nonnegative constant of motion A2 = |z^_p-|-|i^p-|-|z^3p/2. 



B{t) 



V^IjJ — IVa 



(44) 



where z/+ is a solution to the equation (1421) with A = 0. 
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In this case the equations fHT]) . fHHj) and flH^ can be greatly simphfied if 
we put 

u4t) = le\t). (45) 

The result is 



^3 = y(fe2-*ee), 



(46) 



where e{t) satisfies the equation 

e - je + n{t)e = 0, (47) 

m = \fit)\' + Wit) + -2^'H- (48) 

This latter equation admits A2, eq. ( l33l) . as a first integral, 

A2 = ¥ (1 + ^ Ife - + ^ Ife - . (49) 

In f H7|) the term proportional to the first derivative can be eliminated by the 
substitution 

e = e'exp(^i//(r)cir//(r)^, (50) 
which leads to the more simple equation 

e' + Q'{t)e' = 0. (51) 



^'(t) = l/Wr + yit) + 1^ - ^coj + l(^) - f (^). (52) 



It is worth noting at this point that the invariant ladder operators for the 
boson nonstationary forced oscillator have been obtained [9l [10] in terms of 
solutions to the same classical equation (ISTl) . 



5 CS for the fermion forced oscillator 

We define coherent states (CS) for a given fermion system as eigenstates of 
the corresponding invariant fermion annihilation (or creation) operator B{t). 
Since the most general fermion one mode Hamiltonian operator is of the 
form of (nonstationary) forced oscillator ([251), one- mode fermion CS are 
defined as eigenstates of the invariant ladder operator B(t), eq. (H^ (or eqs. 

([26D, my 

Bm;t) = C\C,t). (53) 
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Since B{t) is invariant operator for the FFO, the eigenvalue ( does not depend 
on time t. In terms of the (, B{t), B^{t) and the i?(t)-vacuum |0;t) we have 
for \C;t) the same formulas as for the canonical fermion CS \(), eq. (fTTl) . 
(fT3l) -( |T5|) . in terms of (, b, 6"'" and the 6- vacuum |0). In particular 



\C;t) = e~2^'^{\0;t)-CB\t)\0;t)). (54) 

It remains therefore to construct the (normalized) new ground state |0;t) 
according to its defining equations 

B{t)\0;t) = 0, 

d (55) 
t-\0-t)=Hf\0-t). 

We put 

|0;t) = ao(t)|0) + ai(t)|l), (56) 
substitute this into fl55l) and after some tedious calculations find 



ai{t) = ao{t) ^ 



+ ^ ' (57) 



exp 



where ip^_^ is the phase of i^+(t). The state will represent the exact 
time evolution of an initial canonical CS \() if the initial conditions fl36|) 
are imposed: |C;0) = \(). In this case, as we have shown in section 3, the 
time evolved state \C',t) could be again an eigenstate of b if the oscillator is 
not 'forced', i.e. if f(t) = 0. Let us note that the time-dependence of the 
constructed states is obtained in terms of solutions to the classical system 
(I28l)-( l30|) . or equivalently to the classical equation ( 15T|) . The latter is the 
same equation that appeared in the time evolution of the CS of bosonic FFO 

Our method of construction of dynamical invariants differs slightly from 
the Lewis- Riesenfeld method [12j (developed for bosonic oscillators). Lewis 
and Riesenfeld used to first construct Hermitian invariant, which then is 
represented as a product of normally ordered ladder operators. To make 
connection to their approach let us suppose that we first succeeded to con- 
struct the Hermitian invariant N{t) and to find some ladder operators B{t), 
B\t) that factorize it: N{t) = B\t)B{t). It is clear that 5(t) may differ 
from our non-Hermitian invariant B{t) in a phase factor: B{t) = e^^WS(t). 
We can then in a standard way construct normalized eigenstates of N{t), 

N{t)\0^) = and N {t)\T^) = \T^) , (58) 
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and of Bit), 

m\cvt) = ckTS, (59) 

= (1 - \C0 [m - CIM)] (60) 

which however do not obey the Schrodinger equation since, in general B{t) 
may not be invariant. To obtain solutions |n;t) and the above eigen- 

states |n; t), n = 0, 1, should be multiplied by phase factors, 

|n;t) = e*'^"(*)M), n = 0, 1, (61) 
IC; t) = (1 - \C0 [e'^''^'^ m - Ce^<^^ W | M)] (62) 
which should obey the equations 

= (M^l - HM- (63) 

Evidently the state (l62i) is eigenstate of B{t) with time dependent eigenvalue 
at) = exp(z<^(t)), ^{t) = Mt) - Mi)- 

The phase ip{t) = (pi(t) — 4>Q(t) consists of two parts - geometrical one ip'^, 
and dynamical one (f^ = f — f'^ [17J, 

^""(t) = ^^*((T;V||^|T?f)-('o;7||^|0?f))dt' (64) 
= ip{t) + ^ (J^\H\h^) - {0^\H\0^)^ dt' (65) 

Concluding Remarks 

In this article, we have extend the earlier results of the boson coherence 
Hamiltonian and boson invariant ladder operators to the fermion coherence 
Hamiltonian and fermion invariant ladder operators. We have pointed out 
that unlike the boson coherence Hamiltonian, which is of the more general 
form of the nonstationary forced oscillator, the fermion coherence Hamilto- 
nian is of more simple form of (nonstationary) fermion oscillator. 

For the more general (in fact most general) fermionic system of nonsta- 
tionary forced oscillator we have constructed invariant ladder operators and 
the related Fock and coherent states. We succeeded to express these invari- 
ants and the time evolution of the corresponding states in terms of the same 
classical equation, that describe the evolution of coherent states of the boson 
nonstationary forced oscillator [HI [IH] • The relation of the invariant ladder 
operators method to the Lewis-Riesenfeld method [T2] was briefly described 
on the example of nonstationary fermion systems. 
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